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The notion of entangling power of unitary matrices was introduced by Zanardi, Zalka and Faoro 
[PRA, 62, 030301]. We study the entangling power of permutations, given in terms of a combinatorial 
formula. We show that the permutation matrices with zero entangling power are, up to local 
unitaries, the identity and the swap. We construct the permutations with the minimum nonzero 
entangling power for every dimension. With the use of orthogonal latin squares, we construct the 
permutations with the maximum entangling power for every dimension. Moreover, we show that 
the value obtained is maximum over all unitaries of the same dimension, with possible exception for 
36. Our result enables us to construct generic examples of 4-qudits maximally entangled states for 
all dimensions except for 2 and 6. We numerically classify, according to their entangling power, the 
permutation matrices of dimension 4 and 9, and we give some estimates for higher dimensions. 
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I. INTRODUCTION 

The notion of entangling power of a quantum evolution 
was introduced by Zanardi, Zalka and Faoro [2l| (see also 
HE El). Let H A , H b and H = H A ® H B be Hilbert 
spaces where dim 7^ = dim 7-^ = d. The entangling 
power of a unitary U € U(H) = U(d 2 ) is the average 
amount of entanglement produced by U acting on a given 
(uncorrelated) distribution of product states. As the pure 
entanglement measure we use the linear entropy S'l(-) of 
the reduced density matrix. For \tp) e H, let [IJ 

^l(|^}) -7^r(l-Trp 2 ), where p = T[ B \i>)(rl>\. 
The entangling power of U is defined as [2l| 

e(U) := f f 5 L ([/^ 1 )|V 2 })#1#2, 

(1) 

where dip 1 and dip 2 are normalized probability measures 
on unit spheres. As the linear entropy is a concave func- 
tion of the reduced density matrix, it is an entanglement 
monotone |17j . and therefore a legitimate pure state en- 
tanglement measure. 

Consider now the scenario where Alice and Bob share 
an unknown product state and want to use a unitary 
operator to create a state which is as highly entangled as 
possible. This means that Alice and Bob are looking for 
unitaries with the maximum entangling power. Zanardi 
|22j observed that there are permutation matrices with 
the maximum entangling power over all unitaries in U{TL) 
when d is odd or d = An, but the case d — 4n + 2 was 
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left open. In this paper we extend this result and study 
the entangling power e(P) of a permutation matrix P. 

The paper is organized as follows. In Section II we give 
a combinatorial expression for e(P) (Theorem 2). In Sec- 
tion III we determine the non-entangling permutations 
(Theorem El- In Section IV we construct permutations 
with the maximum entangling power over all unitaries in 
IA{H) for every d ^ 6. Moreover, we construct permuta- 
tions with the minimum entangling power for all d (The- 
orem UJ. As a corollary, we give an upper bound to the 
number of permutations with different entangling power. 
In Section V, we give a complete numerical classification 
of permutation matrices of dimension 4 and 9 according 
to their entangling power. For higher dimensions, we re- 
port numerical estimates. Finally, in Section VI we draw 
conclusions and propose some open problems. 



II. ENTANGLING POWER OF 
PERMUTATIONS 

For our purposes we find it useful to rewrite Equation 
||TJ in another more concrete form as described by Za- 
nardi (see Equation J3J in Lemma (JTJ, below). The 
new expression for e(U) is based on a correspondence 
between quantum operations on bipartite systems and 
certain quantum states which we are now going to recall. 

Let \ip) — J2i j be a state in H and let X be 

the corresponding operator on H, so that X^ = d(i\X\j). 
The singular values of X are equal to the Schmidt coeffi- 
cients of \ip). The state can be written in terms of this 
operator as |V>) = X®I\i/> + ), where — ^7= Yli=i I**) 
is a maximally entangled state and I is the identity oper- 
ator. This relation establishes a bijection between pure 
states in TL = TLa ® T~Lb and operators X acting on 
Ha — Hb- 

Now, following Cirac et al. 4\, suppose Ha and Hb are 
themeselves bipartite: Ha = H\®H2 andHs = Hz®Hi. 
We can repartition H = Ha ® Hb and regard it as the 
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tensor product of Hi ® H3 and W2 <8> W4. Applying the 
above construction to an operator X 13 yields a state 

|*>A|fl = |X)i2|34 != Jf&lS ® ^24)|* + )l3|24, (2) 



where 



l*+> 



13|24 



1 d 

5.?, 



and N is an appropriate normalization factor. This has 
an important physical interpretation: the operator X13, 
acting on the Hilbert space Hai®T~(.B3, corresponds to a 
quantum operation. If this quantum operation is a ten- 
sor product of two quantum operations, one acting on 
the system 1 and another on the system 3, then |-X") 12 |34 
is a product state with respect to the systems A and B. 
In fact, X13 (g) I 2 4 and |-X") i2|34 have the same Schmidt 
decomposition if X13 is a unitary operator |22| . It turns 
out that if both parties A and B share an entangled state 
|A)i2|34 then the non-local quantum operation X13 ®I2A 
can be implemented with certainty by making use of local 
operations and classical communication only (see Equa- 
tion 6 of |4|). With this in mind, the entangling power 
of a unitary U is related to the entanglement of \U) as 
defined in Equation @. We use the linear entropy as an 
entanglement measure. It might clarify matters to give 
the following two extremal examples: 

• (Identity) If U — I then the corresponding state is 
given by 

1 d 

|-0l2|34 = ^ E I") 12 ® l^') 34 ' 

which is separable with respect to the split A\B — 
12 34. In this case we have Sl(\I)) = 0. 

• (Swap) If U = S = Ylij then the corre- 
sponding state is given by 



|S> 



12|34 



1 - 



which is a maximally entangled vector. In this case 
we have S L (\S)) = 1. 

The swap corresponds to the parties interchanging 
their systems, and might therefore be regarded as the 
most non-local operation possible. Nevertheless, the 
swap operation is non-entangling (that is, it never creates 
entanglement). This discrepancy between non-local and 
non- entangling means that we cannot just use the en- 
tanglement of \U) as a measure for the entangling power 
of the matrix U . However with a small modification one 
can still express e(U) in terms of the linear entropy, as is 
done in the following lemma. 



Lemma 1 (Zanardi [22]) The entangling power of 
unitary U G U{TL) is given by 



e(U) 



S L (\U)) + S L (\US))-S L (\S))}, (3) 



d+l 

where < e(U) < ^ @/. 

A permutation of [n] — {1, 2, n} is a bijection from 
[n] to itself. Every permutation p of [n] induces an n x n 
matrix P — (ptj), called a permutation matrix, such that 
Pij = 1 if p(i) = j and pij = otherwise. Equiva- 
lently a permutation on [n] induces a linear map of an 
n-dimensional Hilbert space which permutes a given ba- 
sis of the space (a permutation operator on H.). If n = d 2 
we can replace [n] by [d] x [d] and write p(i, j) = (fey, Uj); 
thus a permutation of [d 2 ] is represented by a pair of d x d 
matrices K = (kij) and L — (hj). The corresponding 
permutation operator permutes the elements \i)\j) of a 
product basis of Ti: 

p{\m) = \h )\h 3 ). 

Theorem 2 Let P = J2ij \^ijhj){ij\ be a permutation 
matrix in U(Tl). The entangling power of P is given by 



e(P) 



Qp — Qps 



d(d-l)(d+l)' 2 



with 



(4) 



where 



aijjn — (lim\ljjn) — djirnj 
biran — \kim\kin/ = &mm* 

The quantity Q ps is the corresponding expression for the 
permutation PS . 

Proof, l From Lemma ^ 



e(P) = 



d* + d 2 -d 4 [Trp 2 PS + Ti-p 2 P } 
d(d- l)(d+ l) 2 ' 



where p P and p PS are the reduced density matrices of 
the states \P)(P\ and \PS)(PS\, respectively. We can 
express G? 4 Tr/5p in terms of the matrices (fey) and (hj). 
Applying the formula in Equation we find the state 
corresponding to P as 

1 d 

|-P)l2|34=-J 2l \ k imi}l2 l S>\limm) 34: . 



This leads to 



p P = Tr B \P)(P\ = ^ E a vm\kmi)(kj m j\ 



3 



By squaring and taking the trace, we obtain 

d 

The same analysis applies to Qps- ■ 

QTPite The problem of classifying bipartite permuta- 
tion operators according to their entangling power now 
reduces to finding what different values can be taken by 
Qp + Qps- Observe that the coefficients <2jj„ and b imn 
(which are either or f ), and the combination 

Tijmn — ^ijmOjijnbimnbjmnj 

which occurs in Q, have the following interpretation: 

• ciijn = 1 if and only if (i, n) and (j, n), which are in 
the same column of the square [d] x [d] , are taken 
by P to elements in the same column; 

• b imn — 1 if and only if (i,m) and (i,n), which are 
in the same row of the square, are taken by P to 
elements in the same row; 

• r ijmn = f if and only if (i,m), (i,n), (j,m) and 
(j, n) , which form a rectangle in [d] x [d] , are taken 
by P to a rectangle with the same orientation. 

We will denote the rectangle (i, to), (i, n), (j, to), (j, n) 
(formed by the intersection of rows i and j and columns 
to and n) as Rij mn ■ We note that if this rectangle is non- 
degenerate, that is i ^ j and to ^ n, then it contributes 
either or 4 to Qp, since 

Tijmn — Tjimn — Tijnrn — Tjinm — Or 1, 

if i = j or m = n but not both, Rijmn contributes or 
2; while if i = j and m = n, then Rijmn contributes 1. 

QTPite 

QTPite 

III. NON-ENTANGLING PERMUTATIONS 

QTPite Two permutation matrices P,Qg U(H) are 
said to be locally unitarily connected (for short, L Li- 
connected) if there are unitaries V acting on Ha and 
W on Hp such that (V <g> W)P = Q. Then V and W are 
actually permutation operators. Note that if two permu- 
tations are LU-connected then they have the same en- 
tangling power. The set of non-entangling permutation 
matrices is denoted by E° . The following result seems to 
be known to specialists (see, e.g., [1 5j). but we provide a 
proof for completeness 

Theorem 3 Let P G IA(H) be a permutation matrix. 
Then P G E if and only if one of the following two 
conditions is satisfied: 



1. P is LU-connected to I; 

2. P is LU-connected to S. 

Proof. The permutation matrix P is non-entangling if 
P(\ip)\x)) is a product state for all \<p) G Ha and |x) G 
Hp- Consider two basis elements and |i)|j2), and 

suppose 

P(\i)\h)) = \i[M), 
P(\i)\h)) = \i'M)- 

Then either i[ = i' 2 or j[ = j' 2 , otherwise P(\i)(\ji) + \j2))) 
would be entangled. Suppose i\ — i' 2 = i! . Then 
for all j3 we must have P(\i)\fz)) — K') Ijs); f° r if n °t 
^(I*)|J3» = |*")|ii> and then P{\i){\j 2 ) + |j 3 ») is entan- 
gled. So there is a permutation pp such that P(\i)\j)) = 
\i')\pp(j)) . Thus as a permutation of the square lattice 
{{hi) '■ 1 < hj < d}, P takes the row («, •) to a row. 
Consider an element (12, ji) in the same column as («, j\). 
The permutation p cannot take («2, ji) to an element in 
the row («',•) , because that is already full. So p must 
take («2,ji) to an element in the column (■, j[). Thus p 
takes the column (•, ji) to a column, ft then follows that 
p takes every row to a row and every column to a column, 
that is P((\i)\j)) — \PA(i))\pp(j))2 for some permutation 
PA- 

In the second case, j[ = j 2 , the permutation p takes 
two elements in a row to elements in a column, and we can 
similarly show that it takes every row to a column and 
every column to a row, that is P(\i)\j)) — \pA{j))\pp{i)) ■ 
■ 

QTPite The probability of sampling a permutation 
matrix in E° over all permutations of [d 2 ] goes to as 
d — > oo. In fact, by Theorem |3| the number of elements 
in E° is 2(dl) 2 , and therefore the probability is 

2(d\) 2 

— =r; > as d — > oo. 



IV. ENTANGLING PERMUTATIONS 

A. Maximum 

QTPite In this section we construct and count the per- 
mutations with the maximum entangling power d/(d+ 1) 
that can be attained by any unitary operator on Ha ® 
Hp. We will make use of latin squares. Recall that a 
latin square of side d is a d x d matrix with entries from 
the set [d] = {1, . . . , d} such that every row and column 
is a permutation of {1, . . . , d}, and two dxd latin squares 
(kij) and (Uj) are orthogonal if (kij,Uj) is a permutation 
of [d] x [d]. Equivalently, two dxd latin squares (kij) 
and (lij) are orthogonal if the ordered pairs (kij,Uj) are 
distinct for all i and j. Euler believed that there are no 
orthogonal latin squares of side 4n + 2. Only in 1960 was 
it shown by Bose, Shrikhande and Parker that, except 
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for side 6, Euler's conjecture was false 0. The study 
of latin squares is an important area of combinatorics 
with connections to design theory, projective geometries, 
graph theory, etc. 0, E3 ■ Orthogonal latin squares 
have been recently applied in quantum information the- 
ory 

Theorem 4 Let P be a permutation operator on Ha <8> 
Hb defined by P(\i)\j)) = Then the entangling 

power of P equals the maximum value e{P) = d/(d+ 1) 
over 14(H) if and only if the matrices {kij) and (Jij) are 
orthogonal latin squares. 

Proof. By Theorem 2, the entangling power of P is max- 
imised when Qp + Qps is minimised. Now Qp is equal 
to the number of rectangles in [d] x [d] which are taken 
to rectangles by P, with the horizontal lines remaining 
horizontal and the vertical lines remaining vertical. This 
is at least dr, since every rectangle consisting of a sin- 
gle point must be taken to a rectangle. It is precisely 
d 2 if and only if no nonzero horizontal line is taken to a 
horizontal line and no nonzero vertical line is taken to a 
vertical line, i.e. if 

kim 7^ ki n whenever m =/= n 

and 

km ^ Ijm whenever i ^ j. 

On the other hand, Qps is equal to the number of rect- 
angles in [d] x [d] which are taken to rectangles by P, 
but with the horizontal lines becoming vertical and the 
vertical lines becoming horizontal. This will be precisely 
d 2 if and only if no nonzero vertical line is taken to a 
horizontal line and vice versa, that is if and only if 

kim 7^ kj m whenever i =/= j 

and 

km 7^ Un whenever m ^ n. 

Together, these are the conditions for the matrices {kij) 
and (Jij) to be latin squares. Since {kij,kj) form a per- 
mutation of {i,j), the two latin squares are orthogonal. 



Corollary For every d ^ 2, 6 there is a permutation ma- 
trix P £ U(H) such that e(P) is maximum over U{H). 

Proof. It follows from Theorem 01 together with the 
fact that there are two orthogonal latin squares for every 
d^2,6 (see, e.g., 0). ■ 



As an example of a permutation matrix satisfying The- 
orem 4, consider 
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For the permutation matrix R we have e{R) = | which 
is the maximum over all unitaries in IA(H) — {Ha <8 Hb) 
where dim Ha = dim Hb = 3. 

By looking at d 2 x d 2 permutation matrices as made 
up of d 2 blocks, we can state an alternative version of 
Theorem 01 

Theorem 5 Let P £ U{H) be a permutation matrix. 
Then e{P) is maximum over U{H) if and only if P sat- 
isfies the following conditions: 

1. Every block contains one and only one nonzero el- 
ement; 

2. All blocks are different; 

3. Nonzero elements in the same block-row are in dif- 
ferent sub-columns; 

4-. Nonzero elements in the same block- column are in 
different sub-rows. 

Proof. By Theorem 2, the quantity Qp is maximum if 
and only if aij m = bi rnn — 1 for all 1 < i,j,n,m < d. In 
this case Qp — d 4 (for example, whenP = /). On the 
other hand Qp is minimum if and only if aij n = Sij and 
bimn = 5mm an d m this case the sum reduces to 

Qp = Z)i,n=i a iinbinn = d 2 (for example, when P = S). 

To obtain the maximum entangling power we need to find 
the permutation matrix that minimizes Qp + Qps- We 
have seen that the maximum entangling power is d/{d + 
1), and in is shown that this value can be obtained 
if and only if 

S L {\U)) = S L {\US)) = S L {\S)) 

or, equivalently, 

Qp = Qps = Qs = d 2 . 

It is then easy to observe that the conditions 1 and 2 
express the minimality of Qp and the conditions 3 and 4 
express the minimality of Qps- ■ 

Theorem[S]implies that in a d 2 x d 2 permutation matrix 
P attaining the maximum value d/{d + 1), every block 
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contains one and only one nonzero entry (as in the above 
permutation matrix R). It is then possible to represent 
P by a d x d array P = (py). The cell Pij specifies the 
coordinates of the nonzero entry in the ij-th block of P. 
For the above permutation matrix P, we have 



R 



11 23 32 
22 31 13 
33 12 21 



Note that the ij-th cell of R is of the form (fcy, kj) where 
K = (kij) and L — (lij) are the orthogonal latin squares 





1 2 3 




1 3 2 


K = 


2 3 1 


and L = 


2 1 3 




3 1 2 




3 2 1 



It follows from Theorem [S] that a permutation matrix P 
has maximal entangling power if and only if P is obtained 
by superimposing two orthogonal latin squares. 
Direct calculations give the following result. 

Theorem 6 The following statements are true: 



For d 
e(P) = 
U{H). 



2 the matrix P = CNOT attains the value 
| which is maximum over all unitaries in 



For d = 6 the value e(P) = §g| is maximum over 
all permutations P G U(TL) and the maximizing P 
is associated to 



11 22 33 44 55 66 
23 14 45 36 61 52 
32 41 64 53 16 25 
46 35 51 62 24 13 
54 63 26 15 42 31 
65 56 12 21 33 44 



(5) 



Proof. Part 1. It has been shown by Rezakhani [l5j that 
for d = 2, the entangling power of any unitary U £11(4) 
is given by 



e(U) 



x {cos(4ci) cos(4c 2 ) 



1 1 

~~ 3 ~ 9 

cos(4ci) cos(4c 3 ) + cos(4c 2 ) cos(4c 3 ))}, (6) 



where c\, Ca,C3 G R and \c%\ < c 2 < c\ < ir/4. It is easy 
to show that e(U) takes its maximum value 4/9 when 
either C\ — c 2 = tt/4, c 3 = or c\ = 7r/4, c 2 = c 3 = 0. 
Then, every permutation matrix which is LU-connected 
to any of the four matrices attains this maximum value 



(see Table I, given in Section V): 
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S-M. 



Part 2. The permutation associated to P arises from 
the latin squares 



12 3 4 5 6 

2 1 4 3 6 5 

3 4 6 5 12 

4 3 5 6 2 1 

5 6 2 14 3 

6 5 12 3 4 



and 



1 2 3 4 5 6 

3 4 5 6 1 2 

2 1 4 3 6 5 
6 5 1 2 4 3 

4 3 6 5 2 1 

5 6 2 13 4 



which are 'as close to being orthogonal as possible' |ll| . 
For P, we have that Qp = 40 and Qps = 36. Since there 
do not exist two orthogonal latin squares of side 6, these 
are the smallest values obtainable, as it is explained by 
the following argument. The array 



11 22 33 44 55 66 
24 13 46 35 62 51 
56 65 12 21 43 34 
63 54 25 16 31 42 
45 36 61 52 14 23 
32 41 53 64 26 15 



represents the action of P on the set [6] x [6]. The ij-th 
cell of P is kl if in P the contribution of the term \kl) (ij\ 
is nonzero (for example, the 22-th cell of P is 13 because 
1 13} (22 1 is nonzero in P). We have printed in boldface 
the symbols that occur twice in a row or in a column. 
Observe that, given any permutation matrix P G [/(36), 
since every element of the set [6] x [6] has to occur exactly 
once in P, it is never possible to have only one pair of 
equal symbols in the same column, without this being 
the case for also another column. Hence, the minimum 
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value attainable by Qp would be 2 steps ahead of 36, 
that is 40. ■ 



We have established a bijection between pairs of or- 
thogonal latin squares of side d and d 2 x d 2 permutations 
with maximal entangling power. The number of pairs of 
orthogonal latin squares of side d is known only for small 
d (see A072377, 0). It is 36 for d = 3 and 3456 for 
d = 4. Note that these values apply to unordered pairs 
of orthogonal latin squares. This means that the number 
of d 2 x d 2 permutations with the maximum entangling 
power is twice the number of pairs of orthogonal latin 
squares of side d. General methods for constructing pairs 
of orthogonal latin squares are presented in |6J, |7j, [14J . 



For a unitary U that reaches the upper bound e(U) = 
d/(d + 1), we have S L (\U)) = S L (\US)) = 1. In other 
words, \U) is maximally entangled with respect to the bi- 
partite splits 12 134 and 23 1 14. It turns out that this is also 
maximally entangled with respect to the split 13 1 24. The 
easiest way to see this is by looking at Equation^defimng 
\U). The state |t/)i3|24 is obtained by two parties A = 13 
and B = 24 sharing a maximally entangled state, by ap- 
plying the unitary U on Alice's system. Since this is a 
local reversible transformation, the entanglement is pre- 
served, hence the state |£/)i 3 | 2 4 is maximally entangled. 
Thus a unitary U is maximally entangling if and only \U) 
is maximally entangled with respect to the three possi- 
ble bipartite splits [2(|. Note that \U) is also maximally 
entangled in all the four splits 1|234, 2| 134, 3|124 and 
4| 123. Then \U) is indeed a maximally entangled state. 
The construction of maximally entangling unitaries thus 
provides us with a canonical way of constructing such 
maximally entangled 4-qudits for all dimensions except 
d = 2 and d = 6. It has been known for some time 
that this is not possible for d — 2, i.e. 4-qubits (see [Io|'). 
leaving only the case d = 6 open. 



B. Minimum 

In this section we construct the permutations with the 
minimum nonzero entangling power that can be attained 
by permutation operators on Ha ® H-b- Two permuta- 
tion matrices P, Q S U(TL) are said to be in the same 
entangling class if e(P) — e(Q). 

Theorem 7 Let P G U(Ti) be a permutation matrix. 
Then e(P) is nonzero but minimum over all permutations 
in U{H) if 



In such a case 
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■ ■■■ Id 
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. dd d(d 


-1) 



e(P) 



8(d-l) 
d{d+l) 2 



(J) 



Proof. Since P is the permutation matrix closest to 
/, it is clear that e(P) is miminum. To find its value, 
we compare Qp with Qj where I denotes the identity 
permutation. The rectangles Rijmn which contribute to 
Qi but not to Qp are those containing (d — 1, d) or (d, d) 
or both, except for the degenerate rectangles contained 
in the bottom line of the square (i.e. those with i = j = 
d). There are (d — l)(d — 2) non-degenerate rectangles 
containing [d— 1, d) but not (d,d — 1), each contributing 
4 to Qi, and d — 1 degenerate vertical rectangles, with 
two equal vertices at (d, d — 1) and two equal vertices at 
(i,d — 1) where i ^ d. Each of these contributes 2 to 
Qi, so the total contribution from rectangles containing 
(d,d - 1) but not {d,d) is 4(d - l)(d - 2) + 2(d - 1). 
There is an equal contribution from rectangles containing 
(d,d) but not (d,d — 1). Finally, there are d — 1 non- 
degenerate rectangles containing both (d, d— 1) and (d, d), 
each contributing 4 to Qj. The total contribution to Qi 
which is not included in Qp is 8(d — l) 2 . Since Qj = d 4 , 
we have 

Q P = d 4 -8(d-l) 2 . 

On the other hand, Qps — d 2 since no nonzero horizontal 
line is taken to a vertical line by P. Hence, by Theorem 
2, e(P) is given by 0. ■ 

Corollary An upper bound to the number of different 
entangling classes of permutations is given by 



B = 2 + -(d 4 - d 2 -8(d-l) 2 ) 



Proof. We can write Qp = 
and to ^ n, the contribution of "the rectangle R 



(8) 



If i ? j 



is 



ij mn 



or 4, which is even. If cither i — j or m — n, the 
contribution of Rij mn is or 2, which is again even. If 

1 = j and m = n, the contribution of Rijmn is 1, and in 
total we have d 2 such terms. It follows that Qp is even if 
and only if d is even. The same analysis applies to Qps- 
Then Qp + Qps is even for all d. 

Now we have seen that the zero entangling power cor- 
responds to Q = d 4 + d 2 , and the maximum to Q = 2d 2 , 
so that \{d 4 — d 2 + 2) is an upper bound to the num- 
ber of classes, where the coefficient 1/2 comes from the 
fact that two consecutives values of Qp (Qps) differ by a 
multiple of 2. We can tighten this bound and obtain the 
value {B — 2) by making use of the fact from Theorem[3 
that the value of Qp + Qps is d 4 - 8(d - l) 2 + d 2 when 
e(P) is nonzero but minimum. The first term 2 on the 
RHS of Equation|S]occurs because of the two classes cor- 
responding to zero and the maximum entangling power. 
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V. NUMERICAL RESULTS 

In this section we report some numerical results. We 
are interested in counting the number of different en- 
tangling classes for different dimensions. We are also 
interested in the average entangling power over all per- 
mutations of a given dimension. The results are given in 
the following tables: 



Entangling Power e(P) 


Number of elements 
in entangling class 



4/9 


8 
16 



TABLE I: Classes of permutations with different entangling 
power and the number of elements in each class for d = 2. 



Entangling Power e(P) 


Number of elements 




in entangling class 





72 


1/3 


2592 


3/8 


864 


5/12 


1296 


182/375 


10368 


23/48 


20736 


1/2 


27432 


25/48 


36288 


13/24 


44064 


9/16 


101376 


7/12 


44712 


29/48 


46656 


5/8 


22464 


2/3 


3888 


3/4 


72 



TABLE II: Classes of permutations with different entangling 
power and the number of elements in each class for d = 3. 



Dimension d 


Number of classes 


Averaf 


;e entangling power 


2 
3 
4 
5 


2 

15 

> 65 

> 190 


4f ~ 0.29 
§ « 0.55 
0.67 ±0.01 
0.74 ±0.01 



TABLE III: Number of classes of permutations with differ- 
ent entangling power and the average entangling power as a 
function of the dimension d. 



VI. CONCLUSION AND OPEN PROBLEMS 

In this paper we have studied the entangling power 
of permutations. We have shown that the permutation 
matrices with zero entangling power are, up to local uni- 
taries, the identity and the swap. For all dimensions, we 
have constructed the permutations with (nonzero) mini- 
mum entangling power. With the use of orthogonal latin 
squares, we have constructed the permutations with the 
maximum entangling power for every dimension. More- 
over, we have shown that this value is maximum over all 
unitaries of the same dimension, with a possible excep- 
tion for 36. Our result enabled us to construct generic 
examples of 4-qudits maximally entangled states for all 
dimensions except for 2 and 6. We have reported numeri- 
cal results about a complete classification of permutation 
matrices of dimension 4 and 9 according to their entan- 
gling power. 



We conclude with a list of open problems: 

• Describe a general classification of the permutation 
matrices according to their entangling power. 

• Give a formula for the average entangling power 
over all permutation matrices of a given dimension. 

• For d = 6, does there exist U G U(H) such that 
e(U) > if? 

• The formula in Equation J3J) only works for the lin- 
ear entropy. It would be desirable to have a similar 
simple formula of the entangling power of a unitary 
in terms of the von Neumann entropy. 

• Study the entangling power of permutation matri- 
ces in relation to multipartite systems. In this con- 
text, it is conceivable that the permutation matri- 
ces with maximum entangling power are related to 
sets of mutually orthogonal latin squares. 
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